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Abstract. We have obtained high quality FORS1/VLT optical spectra of 85 disk H II regions in the nearby spiral 
galaxies NGC 3351, NGC 3521, NGC 4254, NGC 4303, and NGC 4321. Our sample of metal-rich H n regions 
with metallicities close to solar and higher reveal the presence of Wolf-Rayet (WR) stars in 27 objects from the 
blue WR bump (~ 4680 A) and 15 additional candidate WR regions. This provides for the first time a large set 
of metal-rich WR regions. 

Approximately half (14) of the WR regions also show broad C IV A5808 emission attributed to WR stars of 
the WC subtype. The simultaneous detection of C ill A5696 emission in 8 of them allows us to determine an 
average late WC subtype compatible with expectations for high metallicities. Combined with literature data, the 
metallicity trends of WR features and the WC/WN number ratio are discussed. 

The WR regions show quite clear trends between their observed WR features and the H/3 emission line. Detailed 
synthesis models are presented to understand/interpret these observations. In contrast with earlier studies of 
low metallicity WR galaxies, both W(WR) and 7(WR)//(H/3) are here found to be smaller than "standard" 
predictions from appropriate evolutionary synthesis models at corresponding metallicities. Various possibilities 
which could explain this discrepancy are discussed. The most likely solution is found with an improved prescription 
to predict the line emission from WN stars in synthesis models. 

The availability of a fairly large sample of metal-rich WR regions allows us to improve existing estimates of 
the upper mass cut-off of the IMF in a robust way and independently of detailed modeling: from the observed 
maximum H/3 equivalent width of the WR regions we derive a lower limit for M up of 60—90 in the case of a 
Salpeter slope and larger values for steeper IMF slopes. This constitutes a lower limit on M up as all observational 
effects known to affect potentially the H/3 equivalent width can only reduce the observed W (H/3). 
From our direct probe of the massive star content we conclude that there is at present no evidence for systematic 
variations of the upper mass cut-off of the IMF in metal-rich environments, in contrast to some claims based on 
indirect nebular diagnostics. 

Key words. Galaxies: abundances - Galaxies: evolution - Galaxies: ISM - Galaxies: starburst - Galaxies: stellar 
content - Stars: luminosity function, mass function - Stars: Wolf-Rayet 



1. Introduction 

Wolf-Rayet stars (WR) are the descendants of the most 
massive stars. Although they live during a short time 
(Maeder & Conti 1994) these stars have been detected in 
young stellar systems, such as extragalactic HII regions 
(Kunth & Schild 1986) and the so-called WR galaxies 
(Conti 1991, Schaerer et al. 1999b). They are recognized 
by the presence of broad stellar emission lines at optical 
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wavelengths, mainly at 4680 A (known as the blue WR 
bump) and at 5808 A (red WR bump). The blue bump is 
a blend of N v AA4604,4620, N in AA4634,4641, C m/iv 
AA4650,4658 and He n A4686 lines, that are produced in 
WR stars of the nitrogen (WN) and carbon (WC) se- 
quences. In contrast, the red bump is formed only by C IV 
A5808 and it is mainly produced by WC stars. The de- 
tection of these features in the integrated spectrum of a 
stellar system provides a powerful tool to date the onset 
of the burst, and it constitutes the best direct measure of 
the upper end of the initial mass function (IMF). Thus, 
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if WR features are found in the spectra of star forming 
systems, stars more massive than A%r, where AIwr ~ 
25 M for solar metaUicity, must be formed in the burst. 

The IMF is one of the fundamental ingredients for 
studies of stellar populations, which has an important 
bearing on many astrophysical studies ranging from cos- 
mology to the understanding of the local Universe. In par- 
ticular the value of the IMF slope and the upper mass 
cut-off (M up ) strongly influences the mechanical, radia- 
tive, and chemical feedback from massive stars to the ISM 
such as the UV light, the ionizing radiation field, and the 
production of heavy elements. 

A picture of a universal IMF has emerged from numer- 
ous works performed in the last few years (e.g. Gilmore & 
Howell 1998 and references therein). Indeed, these studies 
derive a slope of the IMF close to the Salpeter value for 
a mass range between 5 and 60 Mq. This result seems to 
hold for a variety of objects and metallicities from very 
metal poor up to the solar metallicity, with the possi- 
ble exception of a steeper field IMF (Massey et al. 1995, 
Tremonti et al. 2002). However, the IMF in high metal- 
licity (12+log (O/H) £ (O/H) « 8.92) systems is much 
less well constrained. Different indirect methods to derive 
the slope and M up give contradictory results. 

The detection of strong wind resonance UV lines in 
the integrated spectrum of high metallicity nuclear star- 
bursts clearly indicate the formation of massive stars 
(Leitherer 1998; Schaerer 2000; Gonzalez Delgado 2001). 
In contrast, the analysis of the nebular optical and in- 
frared lines of IR-luminous galaxies and high metallic- 
ity H II regions indicates a softness of the ionizing ra- 
diation field that has beeninterpreted as due to the lack 
of stars more massive than ~ 30 M Q (Goldader et al. 
1997; Bresolin et al. 1999; Thornley et al. 2000; Coziol 
et al. 2001). However, the interpretation of these indi- 
rect probes relies strongly on a combination of models for 
stellar atmospheres and interiors, evolutionary synthesis, 
and photoionisation, each with several potential shortcom- 
ings/difficulties (cf. Garcia- Vargas 1996, Schaerer 2000, 
Stasihska 2002). For example, recently Gonzalez Delgado 
et al. (2002) have shown that the above conclusion could 
be an artifact of the failure of WR stellar atmospheres 
models to correctly predict the ionizing radiation field 
of high metallicity starbursts (see also Castellanos 2001, 
Castellanos et al. 2002b). 

A more direct investigation of the stellar content 
of metal-rich nuclear starbursts has been performed by 
Schaerer et al. (2000, hereafter SGIT00), using the de- 
tection of WR features to constrain M up . They found 
that the observational data are compatible with a Salpeter 
IMF extending to masses M up ij 40 M Q . Most recently, 
a similar conclusion has been obtained by Bresolin & 
Kennicutt (2002, hereafter BK02) from observations of 
high-metallicity HII regions in M83, NGC 3351 and NGC 
6384. 

Here, we present a direct attempt to determine M up 
based on the detection of WR features in metal-rich H n 
regions of a sample of spiral galaxies. To obtain statis- 



tically significant conclusions about M up and the slope 
of the IMF, a large sample of H n regions needs to be 
observed. For coeval star formation with a Salpeter IMF 
and M up =120 M Q at metallicities above solar, ~ 60 to 80 
% (depending on the evolutionary scenario and age of the 
region) of the H n regions are expected to exhibit WR sig- 
natures (Meynet 1995; Schaerer & Vacca 1998, hereafter 
SV98). Thus, to find ^ 40 regions with WR stars (our ini- 
tial aim) a sample of at least 5-7 galaxies with 10 H n 
regions per galaxy needs to be observed. Spectra of high 
S/N (at least 30) in the continuum are also required to 
obtain an accurate measure of the WR features. For this 
propose, we have selected the nearby spiral galaxies NGC 
3351, NGC 3521, NGC 4254, NGC 4303 and NGC4321, 
which have have sufficient number of disk H II regions of 
high-metallicity, as known from earlier studies. 

Our observations have indeed allowed to find a large 
number of metal-rich WR H n regions. The analysis of 
their massive star content is the main aim of the present 
paper. Quite independently of the detailed modeling un- 
dertaken below, our sample combined with additional WR 
regions from Bresolin & Kennicutt (2002) allow us to de- 
rive a fairly robust lower limit on the upper mass cut-off 
of the IMF in these metal-rich environments (see Sect. ||). 

The structure of the paper is as follows: The sample 
selection, observations and data reduction are described 
in Sect. 0. The properties of the H n regions are derived 
in Sect. y. Section [| discusses the trends of the WR pop- 
ulations with metallicity. Detailed comparisons of the ob- 
served WR features with the evolutionary synthesis mod- 
els are presented in Sect. ^[ More model independent con- 
straints on M up are derived in Sect. ||. Our main results 
and conclusions are summarised in Sect. (?]. 

2. Sample selection, observations and reduction 

2.1. Selection of the HII regions 

Our target galaxies (see Table |l|) are selected among 
nearby spiral galaxies where a sufficient number of disk 
H II regions of high metallicity are known from the pre- 
vious studies of Shields ct al. (1991), Oey & Kennicutt 
(1993), and Zaritsky et al. (1994). Inspection of spectra 
from the two latter studies kindly made available to us 
showed that the vast majority of their spectra are not 
deep enough to allow the detection of WR or other stellar 
signatures in the continuum. 

Metallicities 12 + log(0/H) of all known regions were 
estimated from the published [O n] A3727 and [O III] 
AA4959,5007 intensities using the standard R23 "strong 
line" method and various empirical calibrations. For the 
FORS1 multi-object spectroscopic observations described 
below H 11 regions with metallicities above solar (log R23 ^ 
0.6) were given first priority. Secondary criteria taken into 
account in the choice of the known H 11 regions were a 
large H/3 equivalent width, and bright continuum flux at 
~ 4650 A as determined from inspection of the spectra. 
This procedure lead to a first selection of 4 to 7 H 11 regions 



Table 1. Galaxy sample 
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Galaxy 


NED type and activity 


a (J2000) 


<5 (J2000) 


v r 


distance 










[km s" 1 ] 


[Mpc] 


NGC 3351 


SB(r)b, HII Sbrst 


10h43m57.8s 


+Ild42ml4s 


778 


10.0 


NGC 3521 


SAB(rs)bc, LINER 


Ilh05m48.6s 


-00d02m09s 


805 


7.2 


NGC 4254 


SA(s)c 


12hl8m49.5s 


+14d24m59s 


2407 


16. 


NGC 4303 


SAB(rs)bc, HII Sy2 


12h21m54.9s 


+04d28m25s 


1566 


16. 


NGC 4321 


SAB(s)bc, LINER HII 


12h22m54.9s 


+ 15d49m21s 


1571 


15.21 



per galaxy. Other regions with lower metallicities and/or 
lower H/? equivalent widths were retained as secondary 
targets. 

Up to 19 slitlets per exposure can be used for spec- 
troscopy with FORS1. Our primary targets were first po- 
sitioned using the R-band images (see below) and the 
remaining slitlets were filled whenever possible with sec- 
ondary targets. If a slitlet was left without any of our se- 
lected regions, we attempted to target other H II regions 
selected from the Ha images of Hodge & Kennicutt (1983). 
For each galaxy a nuclear spectrum, to be reported upon 
later, was also obtained. 

2.2. Observations 

R band imaging was obtained with FORS1/VLT in april 
2000, and was used to determine the positions of our tar- 
geted regions with sufficient accuracy. Subtracting a local 
average emission from the host galaxy the R band magni- 
tudes of our target H II regions were determined; typical 
magnitudes of rriR ~ 19-21 are found. 

The spectroscopic observations of our sample of H II 
regions were carried out with FORS1/VLT in the second 
2001 trimester. Table || gives informations about the exact 
dates and meteorological conditions during the observa- 
tions. 

The spectral range from 3600 A to 1 /im was covered 
with a "blue" spectrum from 3600 to 6500 A with grism 
300V+10, and a "red" spectrum from 6000 to 10000 A 
with grism 3001+11. The use of a 1" slit width allowed 
to get medium spectral resolution of around 6 A in the 
blue and 12 A in the red. Due to the limited slit size, a 
fraction of the total nebular emission of the regions may 
be lost. This effect is acounted for in our interpretation of 
the data (Sect. ||). Unless WR stars follow systematically 
a different spatial distribution than other stars responsible 
for the continuum emission, a possible loss of continuum 
light does not alter our analysis. 

Exposure times for each galaxy (see Table ||) were 
adapted to obtain in the continuum S/N ~ 30 in the blue, 
(needed for a precise measure of the WR bump) and ~ 
10 in the red (needed to measure the [S in] AA9069,9532 
lines). Spectrophotometric standard stars data were also 
acquired. 



2.3. Data reduction and analysis 

Reduction was carried out using the IRAF and MIDAS 
packages. The first steps consisted in the usual bias sub- 
traction, flatfield division, and 2D wavelength calibration. 
Flux calibration was done using a standard atmospheric 
extinction curve and spectrophotometric standard stars. 
Given that the spectrophotometric standards were not al- 
ways obtained during the night of the observations, we es- 
timate an absolute flux accuracy of ~ 10 %. In addition, 
due to the optimisation for a maximum multiplex, the ob- 
servations were not taken at parallactic angle, leading to a 
slight mismatch between the blue and red spectra. A quan- 
titative analysis of the effetcs of differential refraction has 
not been undertaken here. As the main diagnostics used in 
the present paper lie in a limited wavelength range, and 
the observations have been taken at small airmass, this 
should represent a negligible source of uncertainty. 

For each H II region, a background including sky emis- 
sion and underlying emission from the galaxy was ex- 
tracted from the slitlet sub-image. This procedure was 
non-trivial as this background spectrum had in most cases 
to be determined near the edges of the sub-image, where 
the wavelength calibration may slightly deviate from the 
one of the H n region. Special care has been taken for 
the red spectra, since the sky emission was often several 
times brighter than the H n region emission. We thus re- 
calibrated the background emission spectrum according 
to the H II region by comparing the position (and some- 
times the intensity) of the sky emission lines. This time- 
consuming operation gave very satisfying results and use- 
able spectra up to 1 /im for almost all H n regions. The 
final ID spectra were generally extracted with a 4 " wide 
aperture. 

Line intensities and equivalent width were obtained by 
visually placing a continuum on both sides of the line and 
then integrating all over this range. Errors were estimated 
by moving the continuum upwards by half the value of the 
noise near the line position and re-computing the intensity 
and equivalent width. 

Where possible the following nebular emission lines 
were measured: [O n] A3727, the H Balmer line series in- 
cluding Ha to H9, He i A4471, [O III] A4959,5007, [N n] 
A5201, He i A5876, [O i] A6300, [N n] A6548,6584, He i 
A6678, [S n] AA6717,6731, He i A7065, [Ar in] A7136, [O n] 
A7325, and [S III] AA9069,9532. If present, broad emission 
lines at A ~ 4680 A (referred to subsequently as the (blue) 
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Table 2. Log of the observations with meteorological conditions and exposure times for both grisms 



galaxy 


date 


weather 


seeing ["] 


exp. time blue [s] 


exp. time red [s] 


NGC 3351 


19.04.2001 


photometric 


0.8-1.0 


1700 


1700 


NGC 3521 


25.04.2001 


clear 


1.6-2.0 


1800 


1800 


NGC 4254 


23.05.2001 


clear 


1.1-1.4 


900 


900 


NGC 4303 


23.05.2001 


clear 


0.8-1.1 


750 


750 


NGC 4321 


19.06.2001 


photometric 


1.3-1.5 


1050 


1050 



WR bump), C m A5696, and C iv A5808 indicative of 
Wolf-Rayet (WR) stars were also measured. The spectra 
were also inspected for the presence of stellar absorption 
lines like the Ca n triplet, the CH G band at ~ 4300 A, 
Mg lines at ~ 5200 A, or TiO bands. 

The spectra were deredened using the Whitford et 
al. (1958) extinction law as parametrised by Izotov 
et al. (1994) assuming an underlying absorption of 
W(iip) — 2 A and an intrinsinc Balmer decrement ratio 
of I(Ra)/I{R/3) = 2.86. 

All detailed results including finding charts, line mea- 
surements, and a detailed analysis of the nebular proper- 
ties will be published in a forthcoming paper. 

3. Properties of the HI I region sample 

The properties of our galaxy sample as given by the NED 
database and the adopted distances, are summarised in 
Table |l|. For NGC 3351 and the Virgo cluster member 
NGC 4321 we adopt the Cepheid distances from Freedman 
et al. (2001). The other two Virgo galaxies (NGC 4254, 
NGC 4303) are member of the same subgroup as NGC 
4321 (Boselli, private communication). We therefore adopt 
an identical, approximate distance of 16 Mpc. The dis- 
tance of NGC 3521 is taken from Tully's (1998) Nearby 
Galaxy Catalog. 

A total of 121 spectra were extracted from the 95 
slitlets. Nebular emission lines were detected in 88 spectra; 
85 correspond to extra-nuclear regions. 

3.1. Metallicities 

The metallicity O /H of the H II regions has been estimated 
using the following empirical calibrations: the calibrations 
of Kobulnicky et al. (1999, hereafter KKP) using [O in] 
AA4959,5007/[O n] A3727 and ([O in] AA4959,5007+[O n] 
A3727)/H/3 (—R 23 ) based on the photoionisation model 
grid of McGaugh (1994), the similar P-mcthod of Pilyugin 
(1991), and the older P23 calibrations of Edmunds & Pagel 
(1984) and Zaritsky et al. (1994). 

The O/H abundances obtained from these methods 
are compared in Fig. [l]. Unsurprisingly rather large dif- 
ferences are obtained. As well known, at abundances 
12 + log(0/H) 5s 8.5-8.6 the various R23 methods yield 
similar results, while the differences increase towards 
higher metallicities (see e.g. comparison in Pilyugin 2001). 
Systematically lower values are found from the P-method 
of Pilyugin (2001). Although calibrated only for regions 




12+log(0/H) - KKP calibration 



Fig. 1. Comparison of metallicities O/H of our H II re- 
gions derived from various empirical calibrations. The 
Kobulnicky et al. (1999) calibration is taken as a reference 
(x-axis). Different symbols show O/H derived from the 
Pilyugin P method (filled triangles), and the P 2 3 meth- 
ods of Zaritsky et al. (1994, squares), and Edmunds & 
Pagel (1984, stars). See comments in text. 

with 12 + log(0/H) ^ 8.6, this could indicate a systematic 
overestimate of the absolute metallicities using the other 
methods. To ease comparisons with the recent study of 
BK02 of metal-rich H 11 regions we subsequently adopt 
the KKP calibration by default except otherwise stated. 

The metallicity distribution of our entire sample is 
shown in Fig. |[ The mean metallicity is < 12 + 
log(0/H) >= 8.88 ± 0.22 (8.57 ± 0.24) using the KKP 
(Pilyugin's P) calibration. The vertical dashed line indi- 
cates the solar value (12 + log(0/H)=8.92) adopted in 
McGaugh's calculations used for the calibration of KKP. 

3.1.1. Regions with direct T e determinations 

The transauroral [O 11] A7325 line has been detected in 
11 H 11 regions allowing thus a direct determination of 
the electron temperature from [O 11] A7325/[0 11] A3727. 
Other potential electron temperature indicators, e.g. [S III] 
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Table 3. Statistics of WR regions 



Galaxy 


# blue bump 


# C iv A5808 


# C in A5696 


cand. blue bump 


cand. C iv A5808 


cand. C ill A5696 


NGC 3351 


2 






4 




2 


NGC 3521 


4 


2 


1 


6 


1 


1 


NGC 4254 


9 


8 


1 




1 


1 


NGC 4303 


9 


4 


3 




3 


2 


NGC 4321 


3 




3 


5 


1 


1 


total 


27 


14 


8 


15 


6 


10 




12+log(0/H) - using KKP 13+log(0/H) from KKP 



Fig. 2. Mctallicity distribution of our H II region sam- 
ple (solid line) based on the empirical R23 calibration of 
Kobulnicky et al. (1999). The distribution of O/H for the 
WR regions is shown by the dashed line. The vertical line 
indicates the solar value adopted in the models used by 
these authors. 



Fig. 3. Comparison of metallicities O/H derived from the 
empirical calibration of Kobulnicky et al. (1999, KKP) 
with the determination in regions with measured electron 
temperature T e ([OII]). The errorbars include only the un- 
certainty on the [O 11] A7325/H/3 ratio. The diagonal shows 
the one-to-one relation. Discussion in text. 



A6312, [N 11] A5755, are too weak or could not be mea- 
sured due to the limited spectral reolution. Electron den- 
sities are determined from [S 11] AA6717,6731. T e (0 11) 
and the resulting ionic abundance ratios of ++ /H + and 
+ /H + were derived using this temperature for both ions 
(the atomic data are those listed in Stasinska & Lcithcrcr 
(1996). 

As shown in Fig. || the resulting O/H abundances (as- 
suming O/H = ++ /H + +0 + /H + ) are on average found 
to be lower than those derived from the KKP calibra- 
tion, the largest metallicity being closer to solar. However, 
the O/H derived here are lower limits, due to the strong 
temperature gradients expected at high metallicities (see 
Stasinska 2002). A deeper discussion of the abundances 
in our objects taking into account the observational con- 
straints from the entire emission line spectrum is de- 
ferred to a forthcoming publication. For the purpose of 



the present paper, it is sufficient to note that the bulk of 
our H 11 region sample with low values of R23 have metal- 
licities close to and above solar. 

3.2. Ha luminosities, WR and O star populations 

The histogram of the Ha luminosity of the H 11 regions, as 
measured from our spectra, is shown in the upper panel 
of Fig. U As seen from this figure, ~ 75 % of the H 11 re- 
gions correspond to giant extra-galactic H 11 regions char- 
acterised by L(Ha) ^ 10 38 ergs -1 (Kennicutt 1984, 1991), 
while the remainder are less luminous objects similar to 
normal Galactic H 11 regions. The corresponding number 
of equivalent 07V stars Q A7o7v plotted in the lower panel 
of Fig. [|, ranges from ~ 0.15 07V stars (i.e. presumably 

1 We compute AW = Qo/10 49 ' 12 = L(Ha)/(1.36 10" 12 x 
10 49 ' 12 taking the Lyman continuum flux Qo of an 07V star 
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Fig. 4. Upper panel: Histogram of Ha luminosities of the 
entire H II region sample (solid) and the regions with WR 
detections (dashed). Lower panel: Number of WNL stars 
(derived from blue WR bump assuming WN7 type) as 
a function of the number of equivalent 07V stars de- 
rived from Ha luminosity. Note that the x-axis of both 
plots correspond directly through the definition of N07V 
(see footnote) . Filled (open) symbols indicate regions with 
certain (candidate) WR detections (see text). The lower 
panel shows that secure WR detections are only found for 
regions with -/Vwnl ^ 2-3, as physically expected. This 
indicates that the WR sample in our sample of regions 
with log(Ha) J; 38.2 is fairly complete (cf. text). 



corresponding to ~ 1 late O or early B stars) to ~ 400 
07V stars for the brightest region. 

Our search for WR features in metal-rich H n re- 
gios proved quite successful yielding with 27 WR detec- 
tions a sample of unprecedented size (cf. Castellanos 2001, 
Bresolin & Kennicutt 2002). The number of regions where 
different WR features were detected (hereafter called "WR 
regions" ) at various levels of confidence are listed in Table 
||. The certain WR detections (defined as > 2tr detections) 
are listed in columns 2-4; "candidate" WR regions with 
emission line detections 1.1 < a < 2 are given in cols. 5-7. 
Visual inspection of the spectra yield essentially the same 
detection of the "certain" WR regions. To illustrate the 
quality of our data sample spectra of a secure WR region 
and a candidate region are shown in Fig. ||. 

As also clear from Table [| a large fraction of the WR 
regions shows signatures of WR stars of both WN and 
WC types as anticipated from theoretical expectations 
(Meynet 1995, SV98) and earlier studies of WR galaxies 



Fig. 5. Left panels: FORS1/VLT spectra of a WR region 
in NGC 4254 (top) and a candidate region in NGC 3351 
(bottom) showing also the main line identifications. Right 
panels: Zoom on the spectral region of the blue WR bump 
(top) and the red WR bump (bottom) of the two spectra. 



(Schaerer et al. 1997, 1999ab, Guseva et al. 2000). In our 
sample ~ 50 % of WR regions show WC signatures; pre- 
dictions from the Meynet (1995) and SV98 models yield 
- 30-77 % at metallicities Z ~ 0.008 - 0.040. At least 
1 /3 of the WR regions harbour WC stars of late subtypes 
(WCL), characterised by their strong C in A5696 emis- 
sionPJ The C in A5696/C iv A5808 ratio indicates subtypes 
WC7 or WC8 assuming that the contribution of WN stars 
to C iv A5808 is negligible; if this were not the case the 
mean spectral type could be of later subtype. So far rela- 
tively few WR "galaxies" (true starbursts or extra-galactic 
giant H II regions) with WCL stars are known (cf. Schaerer 
et al. 1999b). However, as late WC types are expected to 
occur preferentially in metal-rich environments (Smith & 
Maeder 1991, Maeder 1991, Philipps & Conti 1992) the 
high detection rate of C in A5696 is not surprising. 

The Ha luminosity distribution of the WR regions is 
shown in Fig. |i| (upper panel, dashed line). Clearly, WR 
stars are only detected in the brightest regions. This is 
not due to the flux limit of our observations as can eas- 
ily be seen by comparison of the smallest WRbump fluxes 
(F(WR) 2ct ~ 4.10~ 16 erg s" 1 cm" 2 ) with the detection 
limit of the faintest emission lines (Fu m (H(3) ~ 10 -16 erg 
s _1 cm -2 ). In fact our observations are essentially deep 



from Vacca et al. (1996) and assuming case B recombination 
for T e = 10000 K and n e = 100 cm" 3 . 



2 It can be seen that C iv A5808 is not formally detected 
in all regions with C in A5696 detections. However, in most 
of these cases (2 exceptions), C iv A5808 is weak, but likely 
present. 
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enough to allow in all galaxies [] the detection of the blue 
WR bump of just ~ 2-3 WNL stars, assuming the average 
4650-4686 A bump luminosity of a WN7 stars of 10 36 5 erg 
s _1 (cf. Smith 1991, Schaerer & Vacca 1998). The num- 
ber of WNL stars derived in this way is plotted in the 
lower panel of Fig. [I| for regions with certain WR detec- 
tions (filled squares) and "candidate" WR regions (open 
circles). 

As our detection limit allows for the detection of few 
(~ 2-3) average WNL stars, the subsample of the bright- 
est regions with F(Hj3) ^ 5.10~ 15 erg s" 1 cm" 2 could 
therefore represent a fairly complete sample of H II regions 
"massive" /bright enough to allow a meaningful compari- 
son between WR detections and non-detections. However, 
a possible bias against regions with small W(H/3) may ex- 
ist (Sect. |). 

In this subsample containing a total of 47 regions we 
find 20 objects without WR signatures, or a fraction of 
~ 57 % regions with WR signatures. Such a high fraction 
of WR detections compares fairly well with the predic- 
tions of 60-80 % by Meynet (1995) and Schaerer & Vacca 
(1998) using the high mass loss stellar evolution tracks at 
metallicities 1/2.5 Z/Zq ^ 2 for bursts with a standard 
Salpeter IMF and an upper mass cut-off M up =120 M . 
Given the fact that very young regions (ages to ~ 1.5-2 
Myr) with large expected H/3 equivalent widths are noto- 
riously absent (in the present sample and other samples 
of H II regions and galaxies) it is, however, not clear how 
significant this finding is. 

4. Trends of WR populations with metallicity 

4.1. Behaviour of the "WR bump" 

Figure |^ shows the WR bump intensities and equivalent 
widths as a function of metallicity for our metal-rich H II 
regions (large filled triangles) and the 11 WR regions in 
spiral galaxies recently reported by Bresolin & Kennicutt 
(2002, large filled squares), together with data compiled by 
Schaerer (1999, small crosses) and Schaerer et al. (2000). 

Our new measurements at high O /H are found to fill in 
the range from the previously observed maximum inten- 
sities/equivalent widths down to lower values. Physically 
the maxima of /(WR) //(H/3) and W(WR) are expected to 
reflect the maximum WR/O star ratio achieved in bursts. 
No lower limit is expected; if present in a given sample, 
such a lower limit presumably reflects the detection limit 
of the WR features. 

The increase of the upper envelope of J(WR)//(H/5) 
with metallicity is known since the work of Arnault et al. 
(1989) and has been reviewed by Schaerer (1999). With 
few exceptions, max(IT(WR)) also seems to show an in- 
crease with O/H as shown here for the first time. The in- 

3 The few remaining "candidate" WR regions at high L(Ha) 
are all located in NGC 4321, have formally a WR detection 
at 1.4 ^cr^2, and appear as rather clear WR detections at 
visual inspection. The observations of NGC 4321 may be of 
less good quality due to limited seeing conditions. 



crease of max(/(WR)//(H / 9)) is naturally interpreted as 
due to the increase of stellar wind mass loss with metallic- 
ity leading to lower minimum mass limit for the formation 
of WR stars, Jl%Ri thereby favouring the presence of WR 
stars at high metallicity (cf. Maeder et al. 1981, Arnault et 
al. 1989, Maeder 1991). Other effects, e.g. a lowering of the 
H/3 flux due to a) increasing amounts of dust absorbing 
ionising radiation or b) lower average stellar temperatures 
at high O/H due to modified stellar evolution, could also 
play a role (cf. Schaerer 1999), but are likely secondary. 

The maxima of the predicted WRbump intensities and 
equivalent widths computed with the code of SV98 with a 
"standard" Salpeter IMF for instantaneous bursts (solid 
line), and extended bursts of duration At = 2 Myr (dot- 
ted), and 4 Myr (long dashed) are overplotted on Fig. [| 
As already shown earlier (cf. Schaerer 1996, 1999, Mas- 
Hesse & Kunth 1999, Guseva et al. 2000) the range of ob- 
servations at subsolar metallicities (12 + log(0/H) ^ 8.6) 
is fairly well reproduced by the models, when account- 
ing for the various uncertainties (e.g. missing H/3 flux in 
slit observations, some objects with small numbers of WR 
stars, some poor spectra; cf. discussion in Guseva et al.). 
The new sample of metal-rich objects plotted here shows 
WRbump strengths smaller than the maxima predicted 
by the "standard" models. The possible reasons for this 
behaviour are discussed in Sect. || where detailed model 
comparisons are undertaken. 

4.2. WC/WN ratio 

We have estimated the relative number ratio of WC and 
WN stars, shown in Fig. [7], in several ways. First the 
number of WN stars, iV(WNL) assuming late WN sub- 
types dominate, is derived from the luminosity of the blue 
WR bump, as described above. The number of WC stars, 
AT(WC), is estimated from the C iv A5808 or C in A5696 
luminosity where measured, again assuming that WN 
stars do not contribute to these lines. As the observed av- 
erage luminosity of WC stars in these lines varies strongly 
with subtype (see SV98), the estimated iV(WC) depends 
on the assumption of the dominant WC subtype. As the 
observations (see above, Guseva et al. 2000, Schaerer et al. 
1999a) indicate that early types (~ WC4) dominate at low 
metallicity, while WC7-8 dominate at high 12 + log(0/H), 
we assume these mean WC subtypes for the sample of 
Guseva et al. (2000). For our high metallicity sample, the 
estimated iV(WC)/iV(WNL) ratios is estimated adopting 
different assumptions on the WC subtype and using C iv 
A5808 or C III A5696 (see Fig. 

The resulting estimates show a fairly clear trend 
of an increasing upper envelope for TV (WC) /TV (WNL) 
with metallicity. Furthermore, and in contrast with the 
limited sample of Guseva et al. (2000), we now find 
at the high metallicity end a number of objects with 
iV(WC)/iV(WNL) > 0.5-1. and a WC/WN number ra- 
tio larger than the observed trend in Local Group galax- 
ies by Massey & Johnson (1998), indicated by the dash- 
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12 + log(0/H) 12 + log(0/H) 

Fig. 6. Observed WR-bump intensities (left panel) and equivalent widths (right panel) as a function of metallicity 
from the compilation of Schaerer (1999, crosses), the samples of Guseva et al. (2000, small filled circles), Castellanos et 
al. (2002a, smalle triangles), SGIT00 (large open squares and circles), BK02 (large filled squares) and the present data 
(large filled triangles). Typical error bars for the Guseva et al. sample are shown. Maximum predicted intensities and 
VF(WR) from the SV98 synthesis models are shown for various star formation histories: instantaneous bursts (solid), 
burst durations At = 2 Myr (dotted) and 4 Myr (long dashed). Note the overall agreement for 12 + log(0/H) 8.5 
and increasing discrepancy above. See discussion in text. 



dotted line in Fig. M. Indeed, while the regions observed by 
these authors are thought to correspond to averages large 
enough to represent the equilibrium iV(WC)/iV(WNL) 
value at constant star formation, larger (and obviously 
also smaller) values should be found in regions with fairly 
short bursts. 

A more quantitative interpretation of the observed WC 
to WN ratio appears difficult for the following reasons. 
First the uncertainties in the estimated iV(WC) /iV(WNL) 
are quite large (cf. above); second, detailed evolutionary 
synthesis model predictions of iV(WC)/iV(WNL) depend 
quite strongly on the adopted interpolation techniques 
(cf. SV98, comparison between results from SV98 mod- 
els and Starburst99 (Leitherer et al. 1999), also Massey 
2002); third, other comparisons with synthesis models re- 
veal potential difficulties (cf. below) . In any case the SV98 
models predict the maximum WC/WN number ratios in- 
dicated in Fig. [j] by the solid line for instantaneous bursts, 
and burst durations of At = 2 Myr (dotted) and 4 Myr 
(dashed) respectively. 

5. Detailed comparison of WR populations with 
synthesis models 

5.1. Procedure 

To interpret quantitatively the observational data we use 
evolutionary synthesis models and proceed essentially as 



in SGIT00. The following main observational constraints 
are used: 

1. Hp and Ha equivalent widths. The former is used as 
a primary age indicator; once W(H(3) is reproduced 
IF(Ha) may serve as an independent consistency test 
for the predicted spectral energy distribution (SED) in 
the red (cf. SGIT00). 

2. Nebular line intensities. F(Ha)/F(H/3) determines the 
extinction of the gas. The use of other line intensities 
requires detailed photoionization modeling which is be- 
yond the scope of this paper. 

3. Intensities and equivalent widths of the main WR fea- 
tures. The blue bump and C IV A5808 (red bump) serve 
as main constraints on the WR population. To avoid 
uncertainties in deblending individual contributions of 
the blue bump we prefer to use measurements for the 
entire bump. In contrast to the spectra of metallic- 
ity objects our spectra show no evident contamination 
from nebular lines (e.g. [Fe in] A 4658, nebular He II). 
To potentially disentangle between various effects (un- 
derlying "non-ionizing" population, loss of photons, 
differential extinction between gas and stars) it is im- 
portant to use both equivalent widths and relative 
7(WR)/7(H/3) intensities (cf. Schaerer et al. 1999a). 

For the model comparisons we use calculations based 
on the evolutionary synthesis code of SV98, which in par- 
ticular includes the most recent calibration of WR line lu- 
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Fig. 7. Estimated number ratio of WC /WNL stars versus 
metallicity. Data derived from C IV A5808/WRbump ob- 
servations from the sample of Guseva et al. (2000) and 
SGIT00 are shown by stars; the assumed mean spec- 
tral type of the WC stars is WC4 for 12 + log(0/H) 
< 8.4 and WC7 for higher metallicities. Results from 
our VLT data are shown for different assumptions: (1) 
from C iv A5808/WRbump assuming a WC7 spectral 
type (filled squares), or WC4 (crosses). (2) from C in 
A5696/WRbump assuming WC7 spectral type (filled 
triangles), or WC8 (open cirles). Maximum predicted 
WC/WN ratios from the SV98 models are shown for in- 
stantaneous bursts (solid), and burst durations At = 2 
Myr (dotted), and 4 Myr (long dashed). The observed 
trend of WC /WN with metallicity in Local Group galax- 
ies, thought to represent an average for constantly star 
forming regions, from Massey & Johnson (1998) is shown 
by the dash-dotted line. 

minosities used to synthesize the WR features, up-to-date 
stellar tracks, CoStar stellar atmospheres for O stars, pure 
H-He models for WR stars and Kurucz models for cooler 
stars (see SV98 for a full description). Except for the im- 
proved O star atmospheres used by SV98 the Starburst99 
synthesis models (Leitherer et al. 1999) use the same basic 
input physics. New generation stellar atmosphere models 
for O and WR stars including a full treatment of non- 
LTE line blanketing and stellar winds have just now be- 
come available for the use in synthesis models (Smith et al. 
2002). However, as the quantities of interest here depend 
only on the total number of Lyman continuum photons 
which is not altered, the use of these more sophisticated 
atmosphere models does not affect our results. 

It is important to stress that in all cases the high-mass 
loss stellar tracks of Meynet et al. (1994) are used. It is 
thought that this adjustment of mass-loss, treated like a 



free parameter, will become ultimately obsolete when a 
proper treatment of the various effects of stellar rotation is 
made in the stellar evolution models. First results tend to 
indicate that this may indeed be the case (Meynet 1999). 
The Meynet et al. (1994) tracks are chosen as they re- 
produce a large number of properties of individual WR 
stars and WR populations (including especially relative 
WR/ O ratios for a standard Salpeter IMF) in Local Group 
galaxies (Maeder & Meynet 1994). The use of other tracks 
(e.g. the "normal" mass loss tracks) which are known to 
disagree with these basic constraints on WR and O star 
populations, would imply a strong inconsistency with the 
Local Group data. 

The basic model parameters we consider are: 

a) Metallicity. Stellar tracks covering metallicities Z = , 
0.008, 0.02 (solar), and 0.04. 

b) IMF slope and upper mass cut-off (M up ). We adopt 
a Salpeter IMF (slope a=2.35), and M up =120 M Q as 
our standard model. 

c) Star formation history (SFH). Models for instanta- 
neous bursts (coeval population), extended burst du- 
rations (constant SF during period At; in this case 
age=0 is defined at the onset of SF, i.e. corresponds to 
that of the oldest stars present), and constant SF are 
considered. 

d) Fraction of ionizing Lyman continuum photons (f-y). 
/ 7 indicates the fraction of ionizing photons absorbed 
by the gas. Our standard value is / 7 = 1. Values / 7 < 1 
are used to simulate various effects (e.g. dust absorp- 
tion, photon leakage outside regions, etc.) leading to a 
reduction of photons available for photoionization. 

Unless stated otherwise our models are calculated as- 
suming an IMF fully sampled over the entire mass range 
(as in SV98). For some cases we have also done model cal- 
culations based on a Monte Carlo sampling of the IMF, in 
order to quantify the effects of statistical fluctuations due 
to the finite number of massive stars. We have verified our 
calculations by comparison with the Monte Carlo models 
and analytical results of Cervino et al. (2000, 2002). 

5.2. Results 

A comparison of the observed equivalent widths and rela- 
tive intensity of the WR bump with standard model pre- 
dictions at different metallicities is presented in Fig. |[ 
The following points can be seen from this Figure: 

— The observed trends of decreasing (increasing) equiv- 
alent width (line intensities) of the WR features with 
decreasing H/3 equivalent width agree with the instan- 
taneous burst model predictions over the same W(Hf3) 
range. Though shorter, part of the initial phase of in- 
creasing (decreasing) FF(WR) (I(WR)//(H/3)) corre- 
sponding to the onset of the WR rich phase does not 
seem to be detected. 

— Most importantly, essentially all the observed WR fea- 
tures are weaker than the predictions of our stan- 




Fig. 8. Observed and predicted equivalent width (left panel) and line intensity with respect to H/3 (right panel) as 
a function of W(Hf3). Our VLT sample is shown by (black) triangles, the BK02 sample with (red) squares. Typical 
uncertainties are 5-10 % for W(H/3), < 10 % for VF(WRbump), and ~ 0.05 dex in log(/(WR)/7"(H/3)). Model predic- 
tions are shown for instantaneous bursts with "standard" IMFs at Z=0.008 (dashed line), Z=Z Q =0.02 (solid line), 
and Z—0.04 (long dashed line). Note the overprediction of the WRbump strength in high metallicity models compared 
with the observations. 



dard models for instantaneous bursts at solar or higher 
metallicity. Very similar observational trends are also 
found in the metal-rich H II region samle of BK02 and 
the metal-rich starbursts of Schaerer et al. (2000). This 
is in stark contrast with observations at lower metal- 
licity where generally an good agreement is found with 
short burst models (e.g. Schaerer et al. 1999a, Guseva 
et al. 2000). 

— The above result is found independently from the ob- 
served WR equivalent widths and relative line intensi- 
ties WRbump/H/3. This indicates that the discrepancy 
between models and observations is not related to the 
possible presence of an underlying older stellar popu- 
lation, which would dilute (reduce) equivalent widths 
but not alter the relative line intensities. 

The following possibilities (one or a combination 
thereof) could be invoked to explain the discrepancy be- 
tween our observations and models: 

1. The metallicities of our HII regions are overestimated. 
Indeed the observations could be reconciled with burst 
models with a "standard" IMF for metallicities Z ~ 
(1/2.5-1) Z©, as shown in the left panel of Fig. || 
(short dashed line). However, despite the uncertain- 
ties in the O/H determinations (cf. Sect. ||) such low 
average metallicities seem very implausible. 

2. Extended bursts. Such a scenario has been invoked by 
SGIT00 for the sample of metal-rich starbursts based 
on the finding of red supergiant features in their spec- 



tra and the fact that these distant objects are mostly 
nuclear starbursts observed through apertures corre- 
sponding to relatively large spatial scales. In this case 
all observed properties could quite well be fitted with 
"standard" solar metallicity models for burst durations 
At ~ 4-10 Myr. However, in view of the different na- 
ture (disk H II regions) of the present sample, indica- 
tions of relatively short formation time scales of H II 
regions (e.g. Massey et al. 1995), and the lack of direct 
signatures of older/red supergiant populations (cf. be- 
low) it seems quite unjustified to appeal to extended 
burst to solve the observed discrepancy. 
3. A modified IMF (upper mass cut-off and/or slope). In 
a plot like Fig. |], a Salpeter IMF with a lower upper 
mass cut-off simply implies that the curve plotted here 
(for M up =12OM ) is joined at lower W(R/3) as the WR 
stars from the most massive stars are absent. This is 
illustrated for the cases of M up = 30 and 60 M Q by the 
shaded domains in Fig. ^. The shape of the predicted 
WR equivalent width or line intensity remains, how- 
ever, unchanged. Therefore the observed discrepancy 
cannot be resolved with an IMF of Salpeter slope and 
a lower value of M up (see also Sect. ||). 
Models with steeper, variable IMF slopes (2.35 < a Si 
3.3) and M up ~ 60-120 M Q could reproduce most of 
the objects, with the exception of the lowest W(H0) 
objects (see Fig. ||). As the least metal-poor objects 
in our sample are probably of similar nature as young 
clusters or H n regions in our Galaxy whose stellar 




Fig. 9. Observed and predicted WRbump equivalent width as a function of W(Hf3). Standard model predictions are 
shown for instantaneous bursts at Z=Zq=0.02 (solid line), and Z=0.04 (long dashed line). Left panel: Thick (green) 
lines with the same styles show models with a standard IMF slope (a — 2.35) and upper mass cut-offs of M up =60 and 
30 Mq delimiting the singly and doubly shaded regions respectively. Right panel: Thick (green) lines with the same 
styles show models with a IMF slope of a = 3.3 and M up =120. 



content has been studied in detail, we may presume 
that their IMF (slope and M up ) should be similar. 
Since none of the Galactic regions have shown convinc- 
ing evidence of a strong deviation of the IMF slope 
from the Salpeter value (see Massey 1998 and refer- 
ences therein), we think that such a steeper slope is an 
unlikely explanation. 
4. Incorrect stellar evolution models and/or "calibration 
data": Although the adopted tracks (Meynet et al. 
1994) compare fairly well with various observations, 
several failures of the non-rotating stellar models are 
also known (see e.g. Maeder 1999). However, the used 
tracks have essentially been calibrated/adjusted to fit 
the observed WR/O ratio in various regions of our 
Galaxy and Local Group objects which are though 
to be at equilibrium, i.e. showing relative popula- 
tions corresponding to constant star formation (see 
compilation in Maeder & Meynet 1994). The relative 
WR/O star ratio is the one most directly related to 
our (time resolved) observables. As this calibration 
yields a fairly good agreement over a large metallic- 
ity range (1/10 5s Z/Zq 2) there seems little room 
for changes in the tracks which could reduce the pre- 
dicted WRbump by the required factor of ~ 2 without 
violating the WR/O constraints in the Local Group. 
One could argue that the calibration data, the ob- 
served WR/O number ratio at solar metallicity and 
above could be incorrect due to possible incomplete- 
ness or biases in the stellar counts (see e.g. related 



discussions in Massey & Johnson 1998). However, to 
reconcile our WR observations in H n regions with the 
corresponding counts for our Galaxy and M31 would 
require a downward revision of the relative WR/O ra- 
tio by up to a factor of 2, which seems highly unlikely. 
5. Uncertainties in synthesis of the WR bump: Presently 
the calculation of observables related to WR stars is 
simply done in the following way in evolutionary syn- 
thesis models. The different emission line strengths are 
computed by multiplying the predicted number of WR 
stars (grouped in different types and/or subtypes) with 
their average line luminosity as derived from observa- 
tions of a sample of WR stars (see SV98). Interestingly 
the intrinsic line luminosity of the strongest line of the 
WR bump, He II A4686, shows a rather large scat- 
ter, namely L4686 = (1-6 ± 1.5) x 10 36 erg s _1 in the 
Galactic and LMC WNL calibration sample of SV98 
with a possible increase of L4686 with the stellar bolo- 
metric luminosity L (see Fig. 1 of SV98). Such a lumi- 
nosity dependence of L4686 with L could in fact (partly 
or fully) explain the observed discrepancy as we will 
now show. 

Splitting the WNL calibration in two domains with 
luminosities above/below logL/L Q = 6, SV98 found 
average line luminosities I/4686 = 5.6 x 10 35 (log L < 6) 
and -L4686 = 3.1 X 10 36 (logL > 6). Replacing in the 
synthesis models the overall average for WNL stars 
by these quantities leads to an important reduction 
of W(WR) in solar metallicity bursts with W(H/3) & 
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Fig. 10. Observed and predicted WRbump equivalent 
width as a function of W(H{3). Standard model predic- 
tions are shown for instantaneous bursts at Z=Zq—0.02 
(solid line), and Z—0.04 (long dashed line). Thick (green) 
lines with the same styles show the predictions with the 
modified L4686 calibration for WNL stars leading to an im- 
portant reduction of the WRbump for logW(H/3) ^ 1.8, 
due to the lower average luminosity of WNL stars in bursts 
with ages Z 4-5 Myr (for Z = 0.02). 



60-70 A, as shown in Fig. 10 4 . At larger H/3 equivalent 
widths (corresponding to ages Ss 4-5 Myr for Z =0.02) 
the WRbump predictions are less modified, since a) 
WC stars contribute more importantly to the bump 
and b) only the youngest bursts with very high W(H/3) 
are dominated by very luminous WNL stars showing 
thus larger W(WR). 

The last option (5) seems the most likely explana- 
tion to explain the surprisingly low WR equivalent widths 
and intensities in our sample of metal-rich H 11 regions. 
Implications on earlier studies of WR galaxies are briefly 
discussed in Sect. 5.2. 



In contrast, the following hypothesis or effects altering 
observed equivalent widths and/or relative line intensities 
cannot be the cause of the discrepancy: 

— Underlying (old) populations diluting the 1F(WR) 
measurements. The line intensities are unaffected by 
underlying populations and show the same discrepancy 
as W(WR) (cf. above). Furthermore, inspection of our 



4 At Z =0.04 one has a reduction of W(WR) both in the 
early burst phase (ages 3.5 Myr) at logW(H/3) <; 1.8 A 
and after ^ 6 Myr (\ogW(R/3) ^ 1.35), while in between the 
models predict that WC stars dominate the bump (cf. SV98, 
Figs. 2 and 11). 



Fig. 11. Equivalent widths of the WRbump as a func- 
tion of the monochromatic continuum luminosity at H/3 
in erg s _1 A -1 (top panel), and as a function of the Ha 
luminosity for the WR regions of our sample. The mean 
and dispersion (1 er) of log L cont (H/3) is plotted in the top 
panel. 



spectra show little or no indication of signatures of 
older stellar populations. 

Differential extinction between gas and stars as fre- 
quently observed in H 11 galaxies and starbursts (cf. 
Fanelli et al. 1988, Calzetti 1997, Mas-Hesse & Kunth 
1999, SGIT00). If present such an effect alters W(RP) 
and 7(WR)//(H/3) but not VF(WR). To bring the ob- 
servations to agreement with standard models would 
require that the gas suffered a lower extinction than 
the stars (implying thus lower corrected VF(H/3) and 
larger 7(WR)//(H/3)), opposite to what is found em- 
pirically. 

Leakage of ionising photons outside the observed re- 
gions, dust absorption of ionising photons, or other 
mechanisms reducing the fraction / 7 of Lyman con- 
tinuum photons used for photoionisation. Correcting 
the observations for such an effect would increase the 
observed W(Hf3) and decrease I(WR)//(H/3), which 
does not reconcile observations and theory (see Fig. 

D- 

Effects due to varying seeing conditions and limited slit 
widths could lead to a loss of nebular emission in the 
aperture or a fraction of the stellar light. The former 
was just discussed ("leakage"). If the WR distribution 
were systematically more extended than the other stars 
responsible of the continuum, this effect could lead to 
reduced WR bump equivalent widths. No such trend 
between W-^(WR) and the seeing conditions is found. 
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Fig. 12. Same as Figs. || (left) and fio| showing the compar- 
ison between model predictions using a fully sampled, an- 
alytical IMF (dotted line, black) and the predicted mean 
FT~(WR) and la scatter (solid line, blue) for a burst of 
scale / mass corresponding to the average observed contin- 
uum luminosity. Observations are shown using the same 
symbols and in the previous Figs. Note, the deviation of 
the mean values for logW(H/3) ~ 1.9-2.2 is due to a nu- 
merical artifact. The comparison shows that no significant 
bias is expected and that the scatter is too small to resolve 
the discrepancy with observations. 



Stochastical fluctuations of the IMF due to small num- 
ber statistics of the massive stars, as invoked by 
Bresolin & Kennicutt (2002). Although indeed ex- 
pected to introduce some scatter (see Cervino et al. 
2000, 2002) this effect cannot explain the discrepancy 
for several reasons. First, the sample discussed here is 
sufficiently large and shows clear observational trends 
with relatively small scatter. In addition, no system- 
atic trend of W(WR) (and I(WR)/J(H/3)) with ab- 
solute scale, such as measured by the continuum lu- 
minosity or Ha luminosity, is observed in our sam- 
ple as shown in Fig. [ll| Second, Monte Carlo simula- 
tions we have undertaken for cluster sizes correspond- 
ing roughly to the observed average continuum lumi- 
nosity of < logL cont >~ 36.8 ± 0.4 f\ predict a typical 
relative scatter of o-(W(WRbump))/W(WR) ~ 25 % 
(cf. Cervino et al. 2002) - too small to account for the 
discrepancy - and no significant bias towards lower 
values as illustrated in Fig. 021 



— The use of different stellar tracks, such as e.g. the 
Geneva tracks with standard mass loss tracks adopted 
by Bresolin & Kennicutt (2002) for the Monte Carlo 
models, which do not reproduce basic constraints of 
massive stars and stellar populations in Local Group 
objects cannot be a solution as they would lead to im- 
portant inconsistencies (cf. discussion in Sect. |5.1| ). 

5.3. Discussion 



In Section 5.2 we have argued that, compared to the nor- 



5 The model L cont at young ages (Js 1 Myr) is scaled to the 
observed value, implying a stellar mass of ~ 5.(3.3) x 10 4 M© 
for a Salpeter IMF between 0.8 (2.) and 120 M . Accounting 
for somewhat older ages would yield a slightly larger mass. 



mal prescription used in our SV98 synthesis models, a 
different prescription should preferrably be adopted to 
predict more accurately the He n A4686 emission from 
WN stars. As several earlier studies including ours (e.g. 
Schaerer 1996, 1999, Schaerer et al. 1999a, Guseva et al. 
2000, SGIT00) are based on the use of the simple average 
He II A4686 line luminosity of SV98 for WNL stars, it is 
important to assess if or to what extent the use of a lu- 
minosity dependent prescription would affect the results 
from previous studies. 

To verify this we have recomputed several sets of 
models for sub-solar metallicities. The maxima of the 
WRbump intensity and VT(WR) (cf. Fig. |^) are only 
slightly modified (increased at 12 + log(0/H) 5s 8.5, and 
decreased above) and lead to a somewhat smaller increase 
with O/H, improving the agreement with the observa- 
tions. For metallicities Z 5s 1/2 Z© the predicted WR 
bump is found to be larger at all ages (as the bulk of WN 
stars are of high luminosity) , whereas for higher metallici- 
ties both larger/smaller WRbump strengths are predicted 
depending on the burst age (W(Hf3)), as for the cases 
shown in Fig. ^| These changes improve the comparison 
with observations at low Z (see e.g. Fig. 7 of Guseva et al. 
2000). No clear statement can be made for intermediate 
metallicities. A better understanding of the dependence of 
the WR emission lines on the stellar parameters appears 
necessary to improve the accuracy of the predictions of 
WR features in evolutionary synthesis models. The im- 
pact of newly available stellar evolution models including 
the effects of rotation on interior mixing and mass loss on 
massive star populations remains also to be explored. 

6. Constraining the upper end of the IMF 

As shown above, and in contrast to previous studies con- 
sidering WR populations in objects with sub-solar metal- 
licities, the quantitative modeling of the WR features in 
metal-rich H n regions is not straightforward. The relia- 
bility of constraints on the IMF obtained from modeling 
of the WR features appears thus unclear at present and 
other constraints are therefore desirable. 

The maximum observed equivalent width of hydrogen 
recombination lines (e.g. W(Hf3)) of a large sample of 
star forming regions provide in principal a constraint on 
the upper mass cut-off of the IMF (e.g Kennicutt 1983, 
Leitherer et al. 1999). In practice, however, as well known 
but not understood to date, there is a notable absence of 



14 



Maximilien Pindao et al.: Metal-rich H II regions and the IMF 




20 40 60 80 100 120 



VI [MJ 

up L 0J 

Fig. 13. Maximum predicted H/3 equivalent width at the 
beginning of the WR phase as a function of M up for so- 
lar metallicity [Z = 0.02)burst models with a Salpeter 
IMF (upper three curves) and a steeper IMF (a = 3.3, 
lower two curves). The dotted curves show models for 
Z = 0.04. The short dashed line corresponds to an ex- 
tended burst of duration At — 4 Myr (Salpeter IMF, 
Z = 0.02.) The observations are plotted at arbitrary M up 
using the same symbols as in Fig. ||. The observed maxi- 
mum (hgWQip) ~ 2.2-2.4) indicates M up ~ 80-90 M Q 
for a Salpeter slope, and ^ 120 M© for a =3.3, or some- 
what lower values for extended bursts. 



regions with W(Hf3) as large as predicted for very young 
bursts with ages between ~ to 1.5-2 Myr. As the onset 
of the WR phase is expected after ~ 2-3 Myr quite inde- 
pendently of the exact adopted stellar tracks (e.g. different 
mass loss scenarios) this sets a new clock, and therefore 
the H/3 equivalent width of the youngest observed WR 
region also contains information on the value of M up . 

In Fig. [l3| we show the dependence of the predicted 
W(Hp) at the beginning of the WR phase (i.e. the max- 
imum of W{HP) during this phase) on the upper mass 
cut-off for different IMF slopes in instantaneous bursts. 
The maximum W(Hf3) depends little on metallicity (see 
the dotted lines) and on the choice of stellar tracks (not 
shown here). Overplotted are the observed VI 7 (H/3) in 
our WR region sample (triangles) and the sample of 
BK02 (squares) drawn at arbitrary M up . The observed 
max(IF(H/3)) (logW/(H/3) ~ 2.2-2.4) indicates an upper 
mass cut-off of ~ 80-90 M Q for a Salpeter IMF or M up ^ 
120 M for a steeper IMF with a =3.3. From all the above 
considerations (Sect.[|) flatter slopes seem excluded. If the 
bulk of the regions were forming stars in extended bursts 
the deduced value of M up has to be lower; for the example 



illustrated here (burst duration At = 4 Myr) this would 
correspond to M up ~ 60 M Q for the Salpeter IMF. 

It is important to note that the value of M up derived 
in this way represents a lower limit. This is the case since 
all observational effects known to affect potentially the 
H/3 equivalent width (loss of photons in slit or leakage, 
dust inside H II regions, differential extinction, underlying 
population) can only reduce the observed W(H{3). The ob- 
served W (H/3) represent therefore lower limits when com- 
pared to evolutionary synthesis models. We thus conclude 
that the available W(H/3) measurements in metal-rich H n 
regions with WR stars yield a lower limit of M up J> 60-90 
Mq for the upper mass-cut off of the IMF. Larger values 
of M up are not excluded. This result is also compatible 
with our favoured models presented in Sect. |E| (see Fig. 
[To|) . Our new estimate of M up , based only on a sample of 
WR regions, provides a more stringent limit than previous 
studies (SGIT00, BK02). 

7. Conclusions 

We have obtained high quality FORS1/VLT optical spec- 
tra of 85 disk H II regions in the nearby spiral galaxies 
NGC 3351, NGC 3521, NGC 4254, NGC 4303, and NGC 
4321. This sample, consisting in particular of a good frac- 
tion of objects with oxygen abundances presumably above 
solar (as estimated from R23 using the calibration reported 
by Kobulnicky et al. 1999), provides an unprecedented op- 
portunity to study stellar populations, nebular properties 
and ISM abundances in H 11 regions at the high metallic- 
ity end. In this hrst paper we have presented the observa- 
tional findings on spectral signatures from massive stars, 
and compared these with evolutionary synthesis models 
with the main aim of constraining the upper part of the 
IMF. 

The average metallicity of our H 11 region sample 
is 12 + log(0/H) ~ 8.9 ± 0.2 using the calibration of 
Kobulnicky et al. (1999). For 12 regions we are able to 
determine the electron temperature from the transauroral 
[O 11] A7325 line, yielding lower limits on O/H (Sect. |). 
For 6 regions we have been able to confirm a high metal- 
licity (12 + log(0/H) ^ 8.8-8.9). Detailed photoionisation 
modeling will be undertaken in the future to improve our 
abundance determinations and to include the full sample 
of H II regions. 

The spectra of a large number (27) of regions show 
clear signatures of the presence of Wolf-Rayet (WR) stars 
as indicated by broad emission in the blue WR bump 
(~ 4680 A). Including previous studies (Castellanos 2001, 
Bresolin & Kennicutt 2002, Castellanos et al. 2002b) 
our observations now nearly quadrupel the number of 
metal-rich H 11 regions where WR stars are known. 
Approximately half (14) of the WR regions also show 
broad C iv A5808 emission attributed to WR stars of the 
WC subtype. The simultaneous detection of C in A5696 
emission in 8 of them allows us to determine an average 
late WC subtype (~ WC7-WC8) compatible with expec- 
tations for high metallicities (Sect. ||). 
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Combined with existing observations of WR regions 
and WR galaxies at sub-solar our data confirm the 
continuation of previously known trends of increasing 
WRbump/H/3 intensity with metallicity, establish also 
such a trend for W(WR), and allow us to estimate the 
trend of the WC/WN ratio with 12 + log(0/H) in extra- 
galactic H II regions (Sect. ||) 

The observed strength of the blue WR bump (relative 
line intensities and equivalent widths) shows quite clear 
trends with W(H./3). Both IV(WR) and 7(WR)//(H/3) 
are found to be smaller than "standard" predictions from 
state-of-the-art evolutionary synthesis models (Schaerer & 
Vacca 1998) at corresponding metallicities. Various possi- 
bilities (including deviations of the IMF from a Salpeter 
slope and a "normal" high upper mass cut-off) which 
could explain this discrepancy have been discussed. The 
most likely solution is found with an improved prescrip- 
tion to predict the line emission from WN stars in syn- 
thesis models (Sect. ||). Using this new prescription the 
observed WR features are found to be broadly consistent 
with short bursts and a "standard" Salpeter IMF extend- 
ing to high masses, as indicated by earlier studies at sub- 
solar metallicities. 

Independently of the difficulties encountered to model 
the WR features in detail, the availability of a fairly 
large sample of metal-rich WR regions allows us to im- 
prove existing estimates (Schaerer et al. 2000, Bresolin & 
Kennicutt 2002) of the upper mass cut-off of the IMF. 
Independently of the exact tracks and metallicity we de- 
rive a lower limit for M up of 60—90 M Q in the case of 
a Salpeter slope, and larger values for steeper IMF slopes, 
from the observed maximum H/3 equivalent width of the 
WR regions. This constitutes a lower limit on M up as all 
observational effects known to affect potentially the H/3 
equivalent width (loss of photons in slit or leakage, dust 
inside H II regions, differential extinction, underlying pop- 
ulation) can only reduce the observed IV (H/3). From our 
probe of the massive star content we therefore conclude 
that there is at present no direct evidence for systematic 
variations of the upper mass cut-off of the IMF in metal- 
rich environments, in contrast to some claims based on 
indirect nebular diagnostics (e.g. Goldader et al. 1997, 
Bresolin et al. 1999, Coziol et al. 2001). What the ori- 
gin of this "universality" of the IMF is, remains an open 
question. 
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